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ABSTRACT 

In this paper, we present a correlation between the spectral index distribution (SED) 
and the dimensionless accretion rate defined as m = iboz/^Bdrf for AGN. This quantity 
is used as a substitute of the physical accretion rate. We select 193 AGN with both 
broad Ha and broad H/3, and with absorption lines near MgIA5175A from SDSS DR4. 
We determine the spectral index and dimensionless accretion rate after correcting for 
both host galaxy contribution and internal reddening effects. A correlation is found 
between the optical spectral index and the dimensionless accretion rate for AGN, 
including low luminosity AGN {Lhoi < lO^'^^erg • sometimes called "dwarf AGN" 
(Ho et al. 1997)). The existence of this correlation provides an independent method 
to estimate the central BH masses for all types of AGN. Wc also find that there is a 
different correlation between the spectral index and the BH masses for normal AGN 
and low luminosity AGN, which is perhaps due to the different accretion modes in 
these two types of nuclei. This in turn may lead to the different correlations between 
BH masses and optical continuum luminosity reported previously (Zhang et al. 2007a), 
which invalidates the application of the empirical relationship found by Kaspi et al. 
(2000, 2005) to low luminosity AGN in order to determine their BLR sizes. 
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1 INTRODUCTION 

Because of the difBculty to calculate the physical accretion 
rate to the BH in an AGN, a dimensionless accretion rate can 
be defined and estimated based on the bolometric luminos- 
ity and BH mass, m = . This parameter can be used 
as the substitute for the actual accretion rate. Moreover, 
the dimensionless accretion rate is an important parame- 
ter in the scheme of the so called unified model for AGN 
(Antonucci 1993, Quintilio & Viegas 1997, Urry & Padovani 
1995), and also the 4D Eigenvector 1 Scheme (e.g., Dultzin- 
Hacyan et al. (2007). The difference in accretion rate leads 
to the principal difference between high luminosity QSOs 
and low luminosity Seyfert galaxies, ft seems to be also one 
of the main physical parameters underlying the 4D Eigen- 
vector 1 Scheme as explained recently in Dultzin-Hacyan et 
al. (2007). In order to obtain the dimensionless accretion 
rate m, two other parameters must be calculated first: BH 
mass Mbh and the bolometric luminosity Lhoi- 

The common method to estimate the bolometric lumi- 
nosity Lboi of AGN is based on the continuum luminosity 
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from the nuclei: Lhoi ~ 9 x L ° given by Kaspi et al. 

5100A & J" 

(2000) and confirmed by Shang et al. (2005) for QSOs. Re- 
cently, the relation was used by Bonning et al. (2007) to 
study the correlation between the accretion disk tempera- 
tures and the continuum colors in QSOs. However, it should 
be stressed that the relation does not hold for ALL kinds of 
AGN, in particular, for low luminosity AGN (Ho et al. 1997a, 
1997b, Ho 1999), because of the different Spectral Energy 
Distribution (SED) (the lack of the Big Blue Bump). 

Several methods are used to estimate BH masses of 
AGN. The most reliable method is based on the stellar veloc- 
ity dispersion of the bulge of the host galaxy first presented 
by Ferrarese & Merritt (2000) and Gebhardt et al. (2000), 
then confirmed by Tremaine et al. (2002) and Merritt & 
Ferrarese (2001) etc. 



Me 



-j^g8.13±0.06 J- 



-)4-»2±»-32Mq (1) 

'200km •s-1'^ " ^ ' 

which indicates a strong correlation between BH masses and 
bulge masses (Haing & Rix 2004, Marconi & Hunt 2003, 
McLure & Dunlop 2002, Laor 2001, Kormendy 2001, Wan- 
del 1999) etc. However we should note that the relation of 
Mbh — cr is obtained through the results of nearby inactive 
galaxies. Whether the relation can be applied to far away 
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active galaxies is an interesting question. So far, there are a 
few dynamical mass estimates of central black holes of broad 
line AGN, and the BH masses axe consistent with the BH 
meisses estimated from the relation of Mbh — <y, although 
the uncertainties are still large. In addition, wc should say 
that the objects in our sample described in the following sec- 
tion are not high luminosity and high redshift QSOs, thus 
the correlation between central black hole and bulge of host 
galaxy can be reasonably considered to hold. 

The other methods arc based on the assumption of viri- 
alization of the Broad Line Emitting Regions (BLRs), or at 
least part of them (Peterson et al. 2004, Onken et al. 2004, 
Sulcntic ct al. 2006, Dultzin-Hacyan et al. 2007). In order 
to calculate the parameter of dimcnsionlcss accretion rate, 
the BH mass is necessary. However, for high luminosity and 
high redshift AGN, it is difficult to measure the stellar ve- 
locity dispersions. The assumption of virialization is applied 
for QSOs. In order to estimate BH masses of QSOs based 
on the assumption of virialization, the most convenient way 
is to use the equation: 



Mbh = f X 



Rblrs X al 



G 



= 2.15 X 10*( 



a"6 



\-!( 5100A nU-BHiit 

3000km •s-1'' 4044erg-s-i'' ® 

(2) 



There are, however, some caveats with this method. First, 
the question of whether the relation Rblrs ~ L*^-^^ » found 

BIOOA 

by Kaspi et al. (2000, 2005) can be applied for all AGN, 
in particular high redshift ones. In an attempt to answer 
this question, we have found that the relation is not valid 
for some special kinds of AGN, such as the low luminos- 
ity AGN (Zhang, Dultzin-hacyan & Wang 2007a, Wang & 
Zhang 2003) and the AGN with double-peaked low ioniza- 
tion emission lines (Zhang, Dultzin-Hacyan & Wang 2007b). 
Second, the estimation of the BH masses of high redshift 
AGN by means of Equation (2) will lead to BH masses larger 
than 1O^°M0 (meaning a > 600km ■ s^^), which leads to 
unreasonable masses of the bulge larger than lO'^^M© (Net- 
zer 2003, Sulcntic ct al. 2006). For this reason, finding an- 
other parameter which can be observationally determined, 
related to the dimensionless accretion rate is an important 
task and is the main objective of this paper. The accre- 
tion rate is determined by two properties: the continuum 
luminosity L_^^ ^ and the BH mass Mbh- The contirmum 
luminosity can'"i)e calculated from the observed spectra £is 
discussed in the the next section. Thus, in order to obtain a 
reliable result, we select Equation (1) to estimate the central 
BH masses of AGN rather than Equation (2). 

The accretion disk model has been widely accepted as 
the standard model for AGN. In the NLTE (Non Local Ther- 
modynamic Equilibrium) accretion disk mode, the generated 
SED (Spectral Energy Distribution) is based on three main 
parameters: BH masses Mbh, accretion rate M and the vis- 
cosity parameter a. An expected result is that there should 
be a correlation between the spectral index and the accre- 
tion rate m. In this paper, we answer the question whether 
the observed spectral index can be used to trace the dimen- 
sionless accretion rate. In section II, we present the data 
sample. Section III gives the results. Finally the discussion 
and conclusions are given in Section IV. In this paper, the 



cosmological parameters Hq = 70km • s ^Mpc ^, f^A = 0.7 
and Q.rn = 0.3 have been adopted. 



2 DATA SAMPLE 

We select objects from SDSS DR4 (Adelman-McCarthy et 
al. 2006) to make up our sample according to the following 
two criteria: First, and most important is that the objects' 
spectra present absorption features, here we focused on the 
absorption line MgIA5175A, in order to measure the stel- 
lar velocity dispersion of the bulge. Second, in order to ob- 
tain the intrinsic continuum luminosity from the rmclci after 
the correction of internal reddening effects using the Balmer 
decrement, Balmer emission lines, both Ha and H/3 must 
be also present. In order to perform accurate measurements 
of the lines mentioned above, several procedures have to be 
followed. 

In order to obtain the continuum luminosity from the 
nuclei we must subtract first the contribution of stellar light. 
An efficient method to subtract the stellar light is the PCA 
(Principle Component analysis) method described by Li et 
al. (2005) and Hao et al. (2005), using the eigenspectra 
from pure absorption galaxies from SDSS or the eigenspectra 
from stars in STELIB (Le Borgne et al. 2003), because the 
method of Principle Component Analysis (PCA) provides a 
better way to constrict more favorable information from a 
series of spectra of stars or galaxies into several eigenspec- 
tra. Here, we used the method from Hao et al. (2005). The 
eigenspectra are calculated by KL (Karhunen-Loeve) trans- 
formation for about 1500 pure absorption galaxies selected 
from SDSS DR4. Then, the first eight eigenspectra and the 
spectra of an A star (which is used to account for star for- 
mation) selected from STELIB (Le Borgne et al. 2003) are 
used to fit the stellar properties of the observed spectra. Af- 
ter this, rather than a power law, a three-order polynomial 
function is used to fit the featureless continuum, because the 
study of composite spectra of AGN shows that the contin- 
uum should be best fitted by two power laws with a break 
of'- 5000A (Francis et al. 1991, Zheng ct al. 1997, Vanden 
Berk et al. 2001). After the last step, the featureless contin- 
uum and the stellar components are obtained based on the 
Levenberg-Marquardt least-squares minimization method. 

After the subtraction of stellar components and the con- 
tinuum emission, the line parameters of emission lines can be 
measured by Levenberg-Marquardt least-squares minimiza- 
tion: one gaussian function for each forbidden emission line, 
two gaussian functions (one broad and one narrow) for each 
permitted emission line. For [OIII] A4959, 5007A, we use an 
extra gaussian function for the extended wings as shown in 
Greene & Ho (2005a). Then we select the objects with re- 
liable broad Ha and broad H/3 according to the following 
criteria: cr{B) ^ 3 x cr{B)err, flux{B) ^ 3 x flux{B)err 
and ij{B) 600km • s~^, where 'B' represents the values for 
the broad Balmer components, 'err' means the measured er- 
ror of the value, cr is the second moment of broad Balmer 
emission lines. 

Then it is necessary to measure the stellar velocity dis- 
persions of the objects selected. However, the accurate mea- 
surement of stellar velocity dispersion is an open question, 
because of the known problems with the template mismatch. 
A commonly used method is to select spectra of several 
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kinds of stars (commonly, G and K) as templates, and then 
broaden the templates by the same velocity to fit stellar fear 
tures, leaving the contributions from different kinds of stars 
as free parameters (Rix & White 1992). However, more in- 
formation about stars included by the templates should load 
to more accurate measurement of stellar velocity dispersion. 
According to the above mentioned method to subtract stel- 
lar components, we created a new template rather than sev- 
eral spectra of G or K stars as templates. Thus, wc apply the 
PGA method for all 255 spectra of different kinds of stars in 
STELIB. Selecting the first several eigenspectra and a three- 
order polynomial function for the background as templates, 
the value of stellar velocity dispersion can be measured by 
the minimum method applied to the absorption features 
around MgIA5175A within the wavelength range from 5100^4 
to 5300l (Zhang, Dultzin-Hacyan & Wang 2007c). Finally, 
we select the objects for which the measured values of stellar 
velocity dispersions are at least three times larger than the 
measured errors. 

Finally, we select 193 AGN with redshift from 0.015 
to 0.25 and with observed featureless continuum luminos- 
ity within the range from 10*^-^^ erg • s"^ to 10'*^ ''^erg • s"^ 
from about 400000 objects classified as galaxies in SDSS 
DR4. The objects have reliable stellar velocity dispersions 
and reliable broad Balmer emission lines. 



3 RESULTS FROM THE DATABASE OF SDSS 

Before proceeding further, a simple discussion about the ori- 
gin of the featureless continuum emission is given. Basically, 
the possibility of nebular emission can be rejected. Accord- 
ing to the luminosity of Recombination lines, the nebular 
continuum emission at 5100A can be simply estimated by 
L ' ~ 0.1 X Lhb, if the electron temperature 

T — 10'* K is accepted. Thus the effects of nebular emis- 
sion can be neglected. Furthermore, we chock the correla- 
tion between continuum luminosity and the luminosity of 
Balmer emission lines found by Greene & Ho (2005b). The 
result is shown in Figure 1. We should note that the con- 
tinuum luminosity and luminosity of Ha are the values bo- 
fore the internal reddening correction as shown in Greene 
& Ho (2005b). The coincident correlation between L 
and Lho ( including the narrow component) indicates that 
our method to subtract the stellar components is reliable 
to some extent. Based on the correlation of Lhcx — L^'^^'^' 

5100A 

for AGN with high continuum luminosity, we can estimate 
the effects of star formation on the continuum luminos- 
ity. For AGN with low luminosity and low redshift, there 
are two components in narrow Ha, one from the AGN 
LhcAGn and the other one from star formation LhcSf 
(Kauffmann et al. 2003). Moreover, we assumed the con- 
tinuum luminosity also includes two components, one from 
the AGN L ' and the other one from star forma- 

5100A,AGN 

tion L - . Furthermore, there is a strong correlation 
between line and continuum luminmosities shown in Figure 



1, Lhoc sf + Lhoc AGN oc (L ' +L ' )*-^". We 

, iioc,OJ- ii<x,Jii^iM V 6100A, SF iWOA,AGN' 

can simply estimate the effects of star formation on contin- 
uum luminosity, if we accept that Lna.SF = s x Lua^AGN 
and Lhoi.agn oc L^'^^'^' (because the relation is bet- 

Z100A,AGN 

ter applied for high luminosity AGN with less effects of star 



formation): 

Lna SF + Lhcy AGN ~ (L ° -\- L ^ )^ 

' «a,^u,iv V 5100A, SF 5100A,AGN' 

Lhu.agn ~ L 



l + s~ (1 + 



5100A,AGN 

L 

5100A.SF si. 157 



(3) 



L 



5100A,AGN 



If we accepted that star-forming regions contribute 65% of 

the narrow Ho flux as described in Kauffmann et al. (2003), 
the parameter of s can be determined as: 

_ Lhci,SF 



LHa,AGN 

0.65 X I/gg.Af 

LHa,B + 0.35 X HHa,N 



(4) 



where 'N' and 'B' represent the narrow component and 
broad component of Ha. The mean value of LHa,B/LHa,N 
is about 3.96 for the objects in our sample. Then we can 
determine that L - IL = ^ 0.12, i.e., the star- 

6100A, SF' SIOOA.AGJV 

forming regions contribute about 10% of the observed con- 
tinuum luminosity. Thus in the following section, we can 
ignore the effects of star fromation. 

In order to obtain a reliable intrinsic continuum shape, 
the internal reddening effects must be corrected. The com- 
mon way to correct them is through the Balmer decrement. 
Here, we assume the intrinsic Balmer decrement as 3.1 for 
Ha and H/3 expected by Gase B recombination (albeit it is 
debatable wether it can be applied to broad lines) with some 
contributions from collisional excitation. Then, the value of 
E(B-V) can be determined from the Balmer decrement: 

Ha 

E(B - V) = -0.97615448 -|- 1.9866313 log(:jj^) (5) 

Hp 

where jg the observed flux ratio. This equation calculated 

from the R-depcndcnt Galactic extinction curve presented 
by Fitzpatrick (1999) can be used to calculate the value of 
E(B-V) simply through the Balmer decrement. In the follow- 
ing of the paper, the continuum luminosity and luminosity 
of Ha are the ones after the correction of BLRs extinction. 

After the correction of internal reddening effects, the 
spectral index can be determined. Here we select three spec- 

F o F o F = 

tral indices: /^w^ and The BH masses are 

5100A 6800A 6800A 

calculated using Equation (1). The internal continuum lumi- 
nosity after the correction of the internal reddening can also 
be calculated. Thus it is easy to check the correlation be- 
tween dimensionless accretion rate and spectral index. Here 
we should notice that the diminsionless accretion rate for 
low luminosity AGN as discussed in the introduction and in 

the next section is also calculated by m = , ~ — , , 

^Edd ^Edd 

although, the bolometric luminosity of low luminosity AGN 
cannot be correctly calculated by Lf,o( ~ 9 x L = as shown 

5100A 

in Ho (1999). However, to some extent, we can accept that, 
if there is also a simple relation Ltoi ^ k x L <• for low 

^ 5100A 
9xL 

luminosity AGN, the calculated — , ^^wa f,^^ )-,g ^ged as a 
substitute of accretion rate, and it is convenient to compare 
the properties of low luminosity and uornial AGN. The cor- 
relations are shown in Figure 2. The Spearman Rank Corre- 
lation Coefficient is 0.68 with P„uii ~ 1-51 x 10"^'', 0.65 with 
Pnuii ~ 1-42 X 10-^'' and 0.67 with P„m ~ 1-94 x 10"^® for 



4 Zhang et al. 



_MooA. ^44ooA g^^^ ^-i4ooA rcspoctivoly. In order to chock the 

6800A 5100A 6800A 

effects of internal reddening, we also show the correlation 
between the dimensionless accretion rate and the spectral 

F 

index j,^^'^'^^ in the right-bottom panel of Figure 2, without 

6800A 

internal reddening correction. The Spearman Rank Correla- 
tion Coefficient for the correlation without reddening correc- 
tion is about 0.56 with P„uii ~ 2.01 x 10~^^. The unweighted 
best fitted results for the correlations between spectral in- 
dices and dimensionless accretion rates axe shown as solid 
lines that correspond to: 



log(- 



log(- 



9 X L 

: 0.56 + 0.18 X log( ^iltloA ) 



= 0.35 -h 0.14 xlog( 



LEdd 

X L 



5100A ■ 



log(- 



6800A 



= 0.92 + 0.33 X log(- 



Lsdd 

9 X L 



(6) 



5100A ■ 



-'Edd 



log( 



F ' 9 X L = 
-^^^){uncorr) = 0.37 + 0.14 x log( -J^mA) 



The last expression is for the correlation between spectral 

F 

index p'^'^'^^^ and dimensionless accretion rate without inter- 

6S00A 

nal reddening correction. Furthermore, we are interested in 
the absolute scatter in the parameter of spectra index which 
can be calculated by: 



N 



(7) 



where Yi, Yija are the measured value of spectra index and 

the fitted value by the equations listed in Equation (6). Fi- 
nally, we can obtain the scatters as follows, A 



0.128, A 



0.124, A 



log( £21004) 

-F6800A 

0.227 and 



iog(i-um4) 

FrjlQQA 



■ 0.123 



±1440(24 J 
-F6800A 



log(J 



uncorr 



-F6800A 

We also show the correlation between the line width of 
broad Ha and the line width of broad H/3 in Figure 3. The 
Spearman Rank Correlation Coefficient is about 0.82 with 
Pnuii ~ 0. The correlation between the line widths of broad 
Balmer emission lines is: 



aHi3B = (1096.56 ± 124.74) x ( 



lO^km • s- 



°i±"-»^km.s-i 
(8) 

where a is the measured value using a gaussian function 
to measure broad emission lines. This correlation is simi- 
lar to the one for QSOs found by Greene & Ho (2005b), 
FWHMng oc FWHMjjl^^° °^. This indicates that the 
measurement of line parameters of broad Balmer emission 
lines is reliable. In addition, it is convenient for us to com- 
pare the two kinds of BH masses estimated from Equation 
(1) and Equation (2). Before proceeding further, we should 
notice that Equation (2) cannot be applied to the low lumi- 
nosity AGN discussed in the next section, because of the un- 
reasonable correlation between the size of BLRs and the con- 
tinuum luminosity found by Kaspi et al. (2000, 2005) (Wang 
& Zhang 2003, Zhang, Dultzin-Hacyan & Wang 2007a). 
Thus, here, we select the normal AGN in our sample to 



estimate the virial BH masses with Equation (2). Finally 
there are 155 objects shown in Figure 4 to compare the two 
kinds of BH masses estimated by Equation (1) and Equa- 
tion (2). The Spearman Rank correlation coefficient is about 
0.36 with P„uH 4.8 x 10~^, after the internal reddening 
corrections. The correlation also indicates that the measured 
stellar velocity dispersions, the measured internal continuum 
luminosities and the line widths are reliable. 



4 DISCUSSION AND CONCLUSIONS 

There are 38 low luminosity AGN with Lna < Iff'^crg ■ s^^ 
(Ho et al. 1997a, 1997b and Ho 1999), which are shown in 
solid circles in Figure 2. From the figure, we can see that 
there is no difference in the correlation between spectral 

9xL 

index and — , giff^ for normal AGN and low luminosity 
AGN. If the bolometric luminosity of low luminosity AGN 

was different from 9 x L ° , all the low luminosity AGN 
would deviate from the correlation for normal AGN, due 
(probably) to different accretion modes. Even if there is a 
different accretion mode for low luminosity AGN (as sug- 
gested by the lack of the big blue bump in the spectra of 
low luminosity AGN as shown in Ho (1999)), the bolomet- 
ric luminosity of low luminosity AGN can also be calculated 
using Lboi ~ fc X igj^QQ^- Otherwise, we could not find the 

same correlation between spectral index and — , ^^^'^^ for 

^ ^Edd 

low luminosity AGN and normal AGN. 

According to the accretion disk model, the output SED 
is the result of the convolution of other parameters as well, 
such as the central BH mass, the viscosity in the disk, and 
the inclination angle. However there is no correlation be- 
tween the spectral index and central BH masses, which is 
shown in Figure 5. The Spearman Rank Correlation Coef- 
ficient is less than 0.1 with Pnuii > 60% for all objects in 
our sample. An interesting result is that there is actually a 
negative trend (anticorrelation) between BH masses and the 
spectral indexes for the 38 low luminosity AGN. The coef- 
ficient is about -0.54 with Pnuii ^ 4.97 x 10"", -0.52 with 
Pnuii ~ 8.62 X lO"'* and -0.55 with P„y,ii ~ 3.98 x 10"* for 
_^iooA^ _44ooA and respectively. 

eSOOA 5100A 6800A 

Because of the positive correlation between the spec- 
tral index and the accretion rate, a negative correlation be- 
tween the spectral index and central BH masses could be 
expected for all AGN. However our results indicate that 
this expectation is only valid for low luminosity AGN. The 
reason is probably related to the correlation between BH 
masses and continuum luminosity. For Normal AGN, there 
is strong correlation between the BH masses and the con- 
tinuum luminosity (Peterson et al. 2004). However for low 
luminosity AGN, this correlation is much weaker (Zhang, 
Dultzin-Hacyan & Wang 2007a). The correlation between 
the central BH masses and the internal continuum luminos- 
ity is shown in Figure 6. The coefficient is about 0.47 with 
Pnull ^ 9-05 X 10"^*^, however, the coefficient is only 0.12 
with Pnull ~ 49% for the 38 low luminosity AGN. The same 
result for low luminosity AGN can be found in Panessa et 
al. (2006). In their paper, they selected all the low luminos- 
ity Seyfert galaxies from Ho, Filippenko & Sargent (1997a, 
1997b), and found that there is NO correlation between 
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the X-ray or optical emission line luminosities (especially 
[OIIIJA5007A line) and BH masses. The result also confirms 
that there is a different accretion mode for normal and low 
luminosity AGN. 

To estimate the effects of the inclination angle of the ac- 
cretion disk is difficult. However under the assumption that 
the narrow line emission region is isotropic, we can check 
the correlation between the continuum luminosity and the 
luminosity of narrow emission lines. If the objects have very 
different inclination angles of the accretion disk, a loose cor- 
relation between the continuum luminosity and the lumi- 
nosity of narrow emission line should be expected. Here we 
show the correlation between L ° and the luminosity of 
narrow Ha in Figure 7. Although it is more common to use 
the [OHI]A5007A as an isotropic estimator of AGN luminos- 
ity, we prefer to use the narrow component of Ha because 
of the following reason. [OHI] emission line frequently can- 
not be fitted by one single gaussian function, because it has 
extended wings (Greene & Ho 2005a). The two components 
are not emitted from the same region. The extended com- 
ponent is probably emitted from the far-side of the BLRs. 
Thus when we fit the [GUI] line, two gaussian functions 
are applied as described in Section H. We thus select the 
narrow component of Ha rather than [OHI] to test the ef- 
fects of inclination angle. A strong correlation can be con- 
firmed. The speaxman Rank Correlation Coefficient is about 
0.89 with Pnuii ~ for normal AGN, and about 0.67 with 
PnuH ~ 4.04 X 10"^*' for the 38 low luminosity AGN. The 
best fit to the correlation (after considering the error in the 
determination of the luminosity of narrow Ha) is given by: 

log(L/f„jv) = (1.373±0.032)-F(0.915±0.003)xlogL ^ ^erg-s"^ 

5100^ (9) 

This result indicates that the effects of the inclination angle 
can be neglected for the correlation between the spectral 
index and the dimensionless accretion rate. 

The correlation between spectral index and dimension- 
less accretion rate found in this research provides another 
independent method to estimate the central BH masses of 
AGN. The spectral index and continuum luminosity can be 
directly determined from the observed spectrum in the op- 
tical band, and subsequently the Eddington Luminosity, i.e. 
BH masses, can be determined by means of the correlation 
we found. This method has the advantage of being inde- 
pendent of the different correlations between the size of the 
BLRs and the continuum luminosity in Equation (2). Also, 
this method can be applied when it is not possible to mea- 
sure the stellar velocity dispersion of the bulge. In future 
work, we will estimate the BH masses of QSOs with higher 
redshift using this method to solve the problem of why virial 
BH masses of QSOs estimated by Equation (2) lead to BH 
masses larger than 10^° Mq, while observational results (for- 
tunately) seem to contradict this result (Dultzin-Hacyan et 
al. 2007). 

Finally, a simple summary is as follows. We first select 

193 AGN with both broad Ha and broad H/3, and with ap- 
parent absorption MgIA5175A from SDSS DR4. Then after 
the determination of the spectral index (after the correc- 
tion of internal reddening effects through the Balmer decre- 
ments for broad Balmer emission lines, and after the sub- 
traction of stellar component) and dimensionless accretion 

rate ( — ^-aiooA )^ find a strong correlation between these 



parameters for AGN, which provides another independent 
and method to estimate the central BH masses of AGN. 
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University of Washington. 
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Figure 1. The correlation between the continuum luminosity and the luminosity of Ha before the internal reddening correction. The 
solid line represents the correlation found by Greene & Ho (2005b). 




Figure 2. The correlation between the spectral index and the dimensionless accretion rate. The solid line represents the unweighted best 
fitted result. Solid circles represent low luminosity AGN with luminosity of Ha less than 10''^ erg • s~^ as identified in Ho et al. (1997a, 
1997b), open circles are normal AGN. 
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Figure 3. The correlation between the line widths of broad Ha and broad H/3. The solid line represents the relation: <THa{B) = angiB). 




Figure 4. The correlation between the two kinds of BH masses, Mbu{(7) estimated from Equation (1) and MguiVirial) estimated 
from Equation (2). The solid line represents the relation; log(AfBj:/ (ct)) = 1.06 log(MB// (^i'**'*'))- Solid circles represent calculations 
after the internal reddening correction, open circles are the values before the internal reddening correction. 
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Figure 5. The correlation between the spectral index and the central BH mass. Solid circles represent the low luminosity AGN with 
luminosity of Ha less than lO^^erg • as identified in Ho et al. (1997a, 1997b). Open circles are the normal AGN. 
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Figure 6. The correlation between the central BH mass and the continuum luminosity after the correction of internal reddening effects. 
Solid circles represent the low luminosity AGN with luminosity of Ha less than lO'^^erg ■ as identified in Ho et al. (1997a, 1997b),open 
circles are the normal AGN. 
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Figure 7. The correlation between the continuum luminosity and the luminosity of narrow Ha after the correction of internal reddening 
effects. Solid circles represent the low luminosity AGN with luminosity of Ha less than 10*^ erg ■ s^^ as identified in Ho et al. (1997a, 
1997b), open circles are the normal AGN. The solid line represents the best fit Lfj^^^ oc i,0-^^5±0.003 ^ 



